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This  is  the  first  of  two  volumes  dealing  with  experimental  and 
theoretical  results  from  the  scattering  of  light  by  deformed  liquid 
droplets  and  droplets  with  inclusions.  Characterizing  droplet 
deformation  could  lead  to  improved  measurement  of  droplet  size  by 
commercial  aerodynamic  particle  sizing  instruments.  Characterizing 
droplets  with  inclusions  can  be  important  in  some  industrial, 
occupational,  and  military  monitoring  scenarios.  The  nozzle 
assembly  from  a  TSI  aerodynamic  particle  sizer  was  used  to  provide 
the  accelerating  flow  conditions  in  which  experimental  data  were 
recorded.  A  helium-neon  (HeNe)  laser  was  employed  to  generate  the 
light-scattering  data,  and  an  externally  triggered,  pulsed  copper 
vapor  laser  provided  illumination  for  a  droplet  imaging  system, 
arranged  orthogonal  to  the  HeNe  scattering  axis.  The  observed 
droplet  deformation  correlates  well  over  a  limited  acceleration 
range,  with  theoretical  predictions  derived  from  analytical  solution 
of  the  Navier-Stokes  equation. 
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1.  Introduction 


The  authors  are  currently  investigating  light-scattering  methods  that  could 
be  used  to  correct  for  the  errors  that  are  known  to  occur  in  the  measurement 
of  liquid  droplets  by  commercial  aerodynamic  particle  sizing  instruments. 
In  the  course  of  the  investigation,  it  became  apparent  that  the  apparatus 
used  for  droplet  deformation  studies  was  also  capable  of  reproducibly 
generating  droplets  with  inclusions.  These  droplets  with  inclusions  are  im¬ 
portant  because  they  can  be  representative  of  aerosols  found  in  hospital  or 
battlefield  environments,  where  biological  organisms  can  survive  for  pro¬ 
longed  periods  when  contained  within  a  protective  liquid  coating.  This  vol¬ 
ume  therefore  describes  the  acquisition  of  experimental  light-scattering  data 
from  both  deformed  droplets  and  droplets  with  inclusions.  These  data  may 
be  used  to  help  develop  and  test  theoretical  inversion  algorithms,  some  of 
which  are  described  in  volume  2,  which  could  ultimately  allow  rapid  char¬ 
acterization  and/or  identification  of  these  droplet  morphologies.  This,  in 
turn,  could  lead  to  improvements  in  commercial  aerodynamic  particle  sizers 
and  new  opportunities  for  the  development  of  light-scattering  aerosol  moni¬ 
toring  instruments  for  bioaerosol  detection. 

1.1  Aerodynamic  Particle  Sizers 

The  size  of  an  airborne  droplet  may  be  defined  in  a  number  of  ways,  de¬ 
pending  on  the  method  of  measurement.  However,  in  studies  of  therapeu¬ 
tic  aerosol  sprays,  combustion  aerosols,  paint  sprays,  and  coatings,  where 
on-line  or  real-time  measurements  are  required  and/ or  where  the  airborne 
behavior  of  the  droplet  is  relevant  to  the  application,  the  aerodynamic  size 
da  is  frequently  sought.  The  parameter  da  is  defined  as  the  diameter  of  a  unit 
density  sphere  with  the  same  gravitational  settling  velocity  as  the  particle 
in  question.  The  traditional  method  of  determining  the  aerodynamic  size  of 
a  particle  was  to  measure  its  gravitational  settling  velocity  with  suitable 
sedimentation  instruments.  This  process  was  slow  and  laborious,  and  com¬ 
mercial  instruments  based  on  the  acceleration  of  particles  were  developed 
as  an  alternative  [1].  The  most  widely  used  commercial  aerodynamic  sizers 
are  the  TSI  aerodynamic  particle  sizer  (APS)  family  of  instruments  (APS 
3300  series,  TSI,  Inc.,  St  Paul,  MN),  which  measure  the  particle  size  by  accel¬ 
erating  a  sample  airstream  in  which  the  particles  are  suspended  and  mea¬ 
suring  the  velocity  acquired  by  individual  particles.  The  larger  particles 
exhibit  greater  inertia  and  thus  accelerate  more  slowly,  thereby  attaining 
lower  velocities.  In  the  TSI  APS  instruments,  an  aerodynamic  accelerating 
nozzle,  in  which  a  sample  flow  is  sheathed  in  filtered  air,  produces  the  nec¬ 
essary  accelerating  flow,  as  illustrated  in  figure  1 .  The  total  airflow  through 
the  nozzle  is  regulated  to  5 1/ min,  4 1/ min  of  which  is  filtered  to  provide  the 
sheath  flow.  The  time  of  flight  of  the  particle  through  two  closely  spaced 
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Figure  1.  Schematic 
diagram  of  particle 
delivery  system  used 
in  TSI  APS  3300  series 
instruments. 
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laser  beams  positioned  a  few  hundred  micrometers  beneath  the  nozzle  exit 
is  recorded  from  the  pulses  of  light  scattered  by  the  particle  to  an  optical 
detector.  The  particle's  aerodynamic  size  is  assessed  by  relating  the  transit 
time  to  an  internal  calibration  function  based  on  measured  times  of  flight  of 
solid  spherical  calibration  particles.  The  size  range  covered  by  this  instru¬ 
ment  is  from  0.5  to  30  pm  aerodynamic  diameter,  depending  on  the  instru¬ 
ment  model  used. 

However,  aerodynamic  sizing  instruments  of  this  type  suffer  a  significant 
limitation  in  that  the  measured  aerodynamic  size  can  be  significantly  af¬ 
fected  by  particle  shape.  When  one  attempts  to  measure  nonspherical  par¬ 
ticles  of  known  density,  the  shape  (and  orientation)  of  each  particle  sub¬ 
jected  to  the  accelerating  air  flow  governs  the  drag  force  it  experiences  and 
hence  affects  the  measured  aerodynamic  size.  The  measurement  of  liquid 
aerosol  droplets  is  subject  to  significant  error  (25%  undersizing  reported  in 
some  cases)  because  the  droplets  deform  to  oblate  spheroids  in  the  acceler¬ 
ating  airflow.  As  a  result  of  this  deformation,  their  cross-sectional  area  in¬ 
creases,  and  they  experience  a  greater  acceleration  than  would  be  the  case 
with  similar-sized  rigid  spheres.  Despite  this  phenomenon  being  well  re¬ 
ported  in  the  past  by  Baron  [2]  and  Griffiths  et  al  [3],  there  is  as  yet  no 
systematic  method  of  measuring  the  degree  of  deformation  experienced  by 
individual  droplets  in  the  instrument-  and  material-specific  calibration 
curves,  derived,  for  example,  through  gravitational  settling  techniques. 
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1.2  Spatial  Light  Scattering 

The  spatial  distribution  of  light  scattered  by  a  particle,  also  referred  to  in 
certain  texts  as  the  two-dimensional  angular  optical  scattering  (TAOS)  pat¬ 
tern  [4],  is  a  complex  function  of  the  size,  shape,  dielectric  structure,  and 
orientation  of  the  particle,  as  well  as  of  the  properties  of  the  illuminating 
radiation  (wavelengths,  polarization  state).  Analysis  of  the  scattering  pat¬ 
tern  can  provide  a  means  of  characterizing  the  shape,  orientation,  and  in¬ 
ternal  structure  of  the  illuminated  particle,  and  many  researchers  have  ex¬ 
ploited  this  property  in  various  ways.  Previous  work  by  the  authors  [5]  has 
explored  the  potential  of  scattering  pattern  analysis  for  particle  shape  and 
size  classification  and  has  demonstrated  [6,7]  how  such  techniques  may  be 
implemented  in  real-time  airborne  particle  measurement  systems  designed 
for  ambient  aerosol  characterization  and  asbestos  fiber  detection. 

Scattering  patterns  can  cover  different  scattering  angle  ranges,  depending 
on  the  light-collection  geometry  used  to  acquire  them.  The  examples  shown 
in  figure  2  (see  figure  caption  for  description  of  particle  types)  were  recorded 
by  the  imaging  of  the  pattern  of  light  scattered  by  individual  particles  onto 
an  intensified  charge-coupled  device  (ICCD)  camera  as  the  particles  were 
carried  by  an  airstream  through  a  laser  beam.  The  particle  transition  time 
through  the  beam  was  ~2  gs,  and  in  each  case  light  scattered  between  5° 
and  30°  scattering  angles  and  throughout  360°  of  azimuth  was  captured  as 
a  256  x  256  pixel  image.  The  beam  direction  is  perpendicular  to  the  paper 
into  the  center  of  each  image.  Each  white  dot  in  the  patterns  corresponds  to 
a  single  scattered  photon,  and  the  images  thus  represent  photon  distribu¬ 
tion  maps  of  several  thousands  to  several  tens  of  thousands  of  scattered 
photons.  The  images  illustrate  the  wide  variations  of  different  particle  shapes 
and  orientations  that  these  patterns  can  assume.  It  was  the  potential  of  spa¬ 
tial  light-scattering  analysis  for  real-time  particle  shape  characterization  that 
initiated  the  fundamental  study  of  droplet  scattering  reported  here.  It  also 
underpins  an  ultimate  aim  of  this  work — to  provide  an  on-line  optical  means 
of  correcting  for  the  errors  in  measured  aerodynamic  size  caused  by  drop¬ 
let  deformation  in  instruments  such  as  those  in  the  APS  3300  series. 
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Figure  2.  Examples  of 
spatial  light  scattering 
patterns  from  indivi¬ 
dual  airborne  part¬ 
icles:  (a)  hematite 
ellipsoid,  2  pm; 

(b)  copper  flake, 

~5  pm;  (c)  sodium 
chloride  crystal; 

(d)  sodium  chloride 
crystal;  (e)  asbestos 
fiber  (crocidolite); 

(f)  irregular  back¬ 
ground  particle; 

(g)  asbestos  fiber 
(chrysotile);  and 

(h)  water  droplet, 

~9  pm  diameter. 
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2.  Apparatus 


An  apparatus  was  constructed  that  allowed  the  acquisition  of  both  spatial 
scattering  patterns  and  real  images  from  individual  droplets  in  the  sub- 
30-pm  size  range  as  they  traversed  the  measurement  space  below  an  APS 
3320  sample  delivery  nozzle  (donated  by  TSI,  Inc.)-  The  apparatus  is  shown 
schematically  in  figure  3.  As  stated  earlier,  the  actual  TSI  APS  3300  series  of 
instruments  incorporates  two  closely  spaced  cross-polarized  beams  and 
measures  particle  time-of-flight  between  the  beams.  However,  for  the  fun¬ 
damental  results  reported  here,  the  spatial  light-scattering  pattern  images 
were  produced  with  a  single  beam  from  a  633-nm,  15-mW  plane  polarized 
helium-neon  (HeNe)  laser  (model  05-LHP-151,  Melles  Griot,  Inc.,  Carlsbad, 
CA).  Furthermore,  while  the  normal  operating  flow  rate  of  the  TSI  APS  in¬ 
strument  is  5  1/min  (1  1/min  of  sample  flow  rate  and  4  1/min  of  filtered 
sheath  flow),  we  have  additionally  observed  droplet  behavior  over  a  range 
of  flow  conditions,  both  less  than  and  greater  than  the  norm,  to  gain  a  greater 
understanding  of  the  morphological  changes  taking  place. 


Figure  3.  Schematic 
diagram  of  apparatus 
used  to  acquire 
scattering  patterns 
and  real  images  from 
individual  airborne 
droplets. 
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2.1  Light-Scattering  Pattern  Acquisition 

The  beam  from  the  HeNe  laser  passes  through  a  quarter-wave  plate  to  ren¬ 
der  the  light  circularly  polarized  before  it  is  focused  through  a  cylindrical 
lens  to  produce  a  beam  of  elliptical  cross  section  (approximately  1.2  mm  by 
100  pm)  at  the  point  of  intersection  with  the  sample  air  flow.  Monodisperse 
droplets  are  generated  into  the  sample  flow  by  a  vibrating  orifice  aerosol 
generator  (VOAG)  (manufactured  by  TSI,  Inc.),  as  described  next.  As  each 
droplet  traverses  the  beam,  light  scattered  in  the  forward  direction  between 
angles  of  5°  and  20°  is  imaged  onto  an  intensified,  asynchronously  trig¬ 
gered  ICCD  camera  (ICCD225,  Photek  Ltd.,  East  Sussex,  UK).  The  lower 
angular  limit  is  set  by  a  beam  stop,  while  the  higher  limit  is  set  to  avoid 
shadowing  of  the  scattered  light  by  the  lower  surface  of  the  sample  deliv¬ 
ery  nozzle.  Images  from  the  camera  are  digitized,  displayed,  and  stored  on 
a  computer  at  a  rate  of  several  images  per  second  for  later  analysis.  A  trig¬ 
ger  signal  for  the  acquisition  of  a  scattering  pattern  by  the  camera  is  de¬ 
rived  from  a  separate  photomultiplier  detector  module,  which  receives  light 
scattered  at  a  higher  scattering  angle  than  the  camera,  as  indicated  in 
figure  3.  The  rising  and  trailing  edges  of  the  signal  from  the  photomulti¬ 
plier  detector  respectively  initiate  and  terminate  the  camera  integration 
period,  which  is  typically  ~2  us  duration. 

2.2  Droplet  Image  Acquisition 

Real  images  of  the  droplets  are  acquired  with  a  charge-coupled  device  (CCD) 
camera  assembly,  mounted  orthogonal  to  the  HeNe  scattering  axis.  In  the 
manner  first  described  by  Baron  [3,8],  the  droplets  are  viewed  against  a 
bright  diffusing  screen.  The  illumination  for  this  screen  is  derived  from  a 
copper  vapor  laser  (Oxford  Lasers  Ltd.,  UK)  that  delivers  38-ns  pulses  at 
wavelengths  of  510  and  578  nm.  Each  pulse  has  an  energy  of  approximately 
2  mj.  Again,  appropriate  triggering  of  the  laser  pulses  is  achieved  by  use  of 
the  signal  from  the  photomultiplier  detector.  The  imaging  system  uses  a 
high  numerical  aperture  microscope  objective  lens  and  a  standard  micro¬ 
scope  tube,  resulting  in  an  overall  image  magnification  of  approximately 
lOx.  A  digital,  nonintensified  CCD  camera  (model  KP-F100,  Hitachi  Denshi 
Ltd.,  Tokyo,  Japan)  was  used  to  acquire  images,  because  it  offered  both  high 
read-out  speed  and  small  pixel  size  (6.7  pm2).  This  pixel  size  produced  im¬ 
ages  of  acceptable  resolution,  the  diameter  of  a  20-pm  droplet,  for  example, 
covering  approximately  30  pixels  within  a  1300  x  1030  pixel  image.  A  Fou¬ 
cault  knife-edge  test  [9]  was  performed  to  determine  the  resolution  of  three 
pixels  for  the  image  quality  attainable  with  this  configuration. 

We  originally  intended  to  acquire  scattering  patterns  and  real  images  si¬ 
multaneously  from  the  same  droplet.  This  would  have  required  appropri¬ 
ate  optical  filtering  to  prevent  the  intense  illumination  from  the  copper  va¬ 
por  laser  from  reaching  the  ICCD  camera.  However,  because  the  pulsed 
radiation  was  broad-field  (as  a  result  of  the  diffusion  screen)  and  multi¬ 
wavelength,  the  task  of  attenuating  this  radiation  to  a  degree  sufficient  to 
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prevent  damage  to  the  ICCD  camera  was  extremely  difficult  without  also 
significantly  attenuating  the  633-nm  scatter  light  flux.  Fortuitously,  the  uni¬ 
formity  of  droplet  size  produced  by  the  VOAG  for  a  given  set  of  generating 
conditions  was  such  that  the  scattering  patterns  and  images  from  succes¬ 
sive  droplets  were  remarkably  reproducible,  the  differences  being  princi¬ 
pally  a  result  of  the  optical  quantization  noise.  We  therefore  decided  to  avoid 
potential  overillumination  of  the  ICCD  by  recording  the  scattering  patterns 
and  real  droplet  images  from  successive  droplets  rather  than  the  same  drop¬ 
let,  ensuring  that  the  ICCD  was  disabled  during  the  firing  of  the  pulsed 
laser. 

2.3  Droplet  Generation  and  Delivery 

The  VOAG  produces  droplets  by  inducing  mechanical  instability  in  a  fine 
liquid  jet  through  excitation  with  a  piezoelectric  oscillator.  When  first  formed, 
droplets  are  comparatively  large  and  compose  a  solute — in  our  case,  either 
oleic  acid  or  glycerol — in  solution  with  AR-grade  isopropyl  alcohol.  Dry 
(<1  %RH)  compressed  air  is  used  to  transport  the  droplet  aerosol  from  the 
VOAG  drying  column  to  the  light-scattering  instrument,  during  which  time 
the  alcohol  evaporates,  to  leave  pure  oleic  acid  or  glycerol  droplets.  A  bleed 
valve  controls  the  total  flow  rate  through  the  nozzle  subsystem.  For  the 
experiments  described,  the  total  flow  rate  was  set  before  the  scattering  data 
were  acquired,  and  the  sheath  air  flow  rate  was  measured  before  and  dur¬ 
ing  acquisition  with  a  bubble  flowmeter.  The  measured  sample  flow  rate, 
consistently  20  percent  of  the  total  flow  rate,  was  found  to  be  stable  to  within 
1  percent  during  the  period  of  experimental  data  acquisition. 
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3.  Experimental  Data 


Figure  4.  Spatial  light¬ 
scattering  patterns 
from  individual 
droplets  as  function  of 
droplet  size  and 
sample  flow  rate. 


All  experimental  data  presented  here  were  recorded  from  oleic  acid  drop¬ 
lets  (with  the  exception  of  those  shown  later  in  fig.  8b,  which  shows  scatter¬ 
ing  from  a  glycerin  droplet).  Figure  4  illustrates  the  changes  that  occur  in 
the  spatial  light-scattering  patterns  from  individual  droplets,  as  a  function 
of  both  droplet  size  and  increasing  sample  flow  rate.  Each  scattering  pat¬ 
tern  was  recorded  over  a  2-gs  period  as  the  droplet  traversed  the  HeNe 
beam  beneath  the  aerosol  delivery  nozzle.  Each  scattering  pattern  therefore 
represents  the  forward  scattering  from  a  deformed  droplet  illuminated  at 
90°  to  its  axis  of  symmetry.  The  images  on  the  top  row  were  recorded  from 
8-pm-diameter  droplets  at  total  flow  rates  from  1  to  6  1/ min.  Subsequent 
rows  show  the  scattering  behavior  of  12-,  16-,  20-,  25-,  and  30-|o.m-diameter 
droplets  for  similar  values  of  flow  rate.  The  gain  of  the  image  intensifier 
was  reduced  for  larger  droplet  sizes  (and  lower  droplet  velocities)  to  mini¬ 
mize  optical  saturation  effects,  although  these  are  still  visible  as  extended 
solid  white  areas  in  some  of  the  larger  droplet/lower  velocity  patterns. 
Ananth  and  Wilson  [10]  have  carried  out  a  theoretical  analysis  of  the  gas 
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and  particle  velocities  present  in  the  APS  delivery  nozzle.  Their  results  pre¬ 
dict  that  at  the  normal  APS  operational  flow  rate  of  5 1/  min,  the  gas  velocity 
at  the  measurement  point  is  approximately  160  m/ s,  and  the  velocity  of  a 
10-|im  spherical  particle  is  approximately  50  m/s. 

At  the  lowest  flow  rate  and  smallest  droplet  size  (top  left  image  in  fig.  4), 
the  scattering  pattern  correlates  closely  with  that  predicted  by  the  Mie  theory 
[11]  for  a  perfect  sphere.  No  droplet  deformation  is  discernible  from  the 
light-scattering  data  in  this  image.  For  increasing  flow  rate  and/ or  droplet 
size,  however,  the  degree  of  droplet  deformation  would  be  evident  from  an 
increasing  ellipticity  of  the  scattering  maxima  and  minima  evident  on  the 
patterns.  This  behavior  appears  well  behaved  to  a  point,  such  as  at  16  pm 
x  5 1/min,  when  bright  horizontal  barring  occurs  on  the  image.  This  is  pre¬ 
dicted  by  the  theoretical  analysis  described  in  volume  2  of  this  report  [12]. 
This  barring  may  occur  either  as  a  predominant  central  white  bar  or  as  a 
series  of  parallel  horizontal  scattering  fringes,  such  as  those  visible  at  20  pm 
x  5  1/min.  Increasing  the  flow  rate  or  droplet  size  further  results  in  more 
complex  interference  patterns  within  the  images,  such  as  those  seen  at  30  pm 
x  5 1/ min.  This  behavior  appears  to  indicate  that  the  droplets  are  undergo¬ 
ing  a  transition  from  simple  oblate  spheroidal  form  to  the  more  complex 
morphologies  first  observed  by  Baron  [2].  It  also  indicates  that  the  task  of 
attempting  to  rapidly  assess  the  deformation  of  a  single  droplet  by  analysis 
of  its  light-scattering  pattern  (as  we  desire  for  on-line  correction  of  defor¬ 
mation  effects  in  aerodynamic  sizing)  is  nontrivial  for  droplets  deformed 
beyond  a  certain  degree. 

Figure  5  shows  light-scattering  pattern  images  for  undistorted  droplet  di¬ 
ameters  of  20,  25,  and  30  pm,  and  total  flow  rates  from  1  to  6  1/ min.  Adja¬ 
cent  to  each  scattering  pattern  is  the  corresponding  droplet  image.  These 
latter  images  indicate  the  changes  in  droplet  shape  that  produce  the  increas¬ 
ingly  complex  scattering  data.  The  essentially  perfect  spheres  undergo  flat¬ 
tening  into  oblate  spheroids,  and  then  become  biconcave  disks  before  form¬ 
ing  menisci  with  concave  upper  surfaces.  This  morphological  transition  co¬ 
incides  with  the  dominant  feature  of  the  respective  scattering  patterns, 
changing  from  a  horizontal  to  a  vertical  bright  band,  the  former  possibly 
the  result  of  a  morphology-dependent  resonance  in  the  droplet,  and  the 
latter  a  function  of  the  increasingly  narrow  horizontal  cross  section  that  the 
droplet  presents.  We  are  currently  investigating  these  effects  further,  with 
the  hope  of  elucidating  the  causes  of  the  dominant  scattering  features. 
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Figure  5.  Spatial  light¬ 
scattering  patterns 
and  droplet  images 
for  20-,  25-,  and 
30-pm-diameter 
droplets  at  various 
sample  flow  rates 
through  APS  aerosol 
delivery  nozzle. 
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Direction  of  patricle  motion 


4.  Quantifying  Droplet  Deformation 


Inversion  of  the  light-scattering  data  to  yield  droplet  deformation  descrip¬ 
tors  is  addressed  in  volume  2  of  this  report  [12].  As  a  comparative  measure, 
a  first-order  estimate  of  droplet  deformation  under  given  flow  conditions 
may  be  obtained  from  the  Feret  ratio  of  the  droplet  image,  defined  as  the 
ratio  of  Feret  diameters,  i.e.,  the  maximum  distances  between  particle  edges 
in  the  x  and  y  directions  of  the  image  (orthogonal  and  parallel  to  the  par¬ 
ticle  trajectory,  respectively).  The  edges  of  the  droplets,  E(x,y),  were  deter¬ 
mined  by  a  computer-based  Gaussian  edge-detection  algorithm  [13].  The 
threshold  value  for  this  algorithm  can  be  set  manually  or  left  to  be  deter¬ 
mined  automatically  from  the  background  intensity.  Since  the  background 
intensity  is  different  for  each  image,  we  used  the  automatic  threshold  func¬ 
tion.  In  addition,  the  graininess  of  the  background  can  make  the  algorithm 
detect  erroneous  edges  there.  However,  the  algorithm  generally  provides  a 
continuous  line  that  is  evident  as  an  edge  in  the  binary  file  that  the  algo¬ 
rithm  creates.  The  Feret  ratio  can  then  be  determined  from  this  binary  file. 
Figure  6  shows  a  plot  of  Feret  ratio  versus  flow  rate  for  20-pm  oleic  acid 
droplets.  As  can  be  seen  from  the  figure,  the  Feret  ratio  is  well  behaved  up 
to  a  4.5-1  / min  flow  rate.  The  reduction  in  the  rate  of  increase  of  the  Feret 
ratio  immediately  after  this  point  coincides  with  the  droplet  beginning  to 
assume  a  meniscus  form  (see  fig.  5),  for  which  the  Feret  ratio  is  less.  The 
change  in  the  curve  at  the  4.5-1/min  flow  rate  coincides  with  the  point  at 
which  the  horizontal  banding  appears  in  the  scattering  pattern.  A  signifi¬ 
cant  part  of  the  experimental  uncertainty  in  the  Feret  ratio  measurements  is 
caused  by  a  blurring  of  the  droplet  images  in  the  direction  of  motion,  espe¬ 
cially  at  higher  flow  rates.  Methods  to  minimize  motion  blur  by  reduction 
of  the  laser  pulsewidth  are  currently  being  explored. 

Liquid  droplet  distortion  produced  in  an  accelerating  flow  field  has  been 
computed  through  analytical  solution  of  the  Navier-Stokes  equation  by 
Bartley  et  al  [14],  and  this  has  led  to  a  means  of  correlating  observed  droplet 
deformation  to  such  parameters  as  droplet  diameter,  viscosity,  surface  ten¬ 
sion,  and  density.  A  high-Reynolds-number  empirical  approximation  to  the 
pressure  external  to  the  droplet  has  been  used  in  these  calculations.  This 
procedure,  described  in  detail  by  Bartley  et  al,  correctly  predicts  the  change 
in  droplet  shape  with  increasing  flow  rate  from  spherical  to  spheroidal  to 
biconcave  disk.  For  example,  figure  7  shows  the  image  of  a  20-pm  nominal 
diameter  oleic  acid  droplet,  together  with  the  corresponding  computed  cross 
section.  Beyond  this  level  of  distortion,  the  theoretically  computed  shape 
pattern  is  not  reliable. 
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Figure  7.  Image  of 
20-jam  nominal 
diameter  oleic  acid 
droplet,  together  with 
corresponding 
computed  cross- 
sectional  shape. 
Computation  used 
density  =  0.90  g/cm3, 
viscosity  =  0.256 
Poise,  and  surface 
tension  =  32  mj/m2. 
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5.  Droplets  with  Inclusions 


During  the  process  of  collecting  experimental  light-scattering  and  particle 
image  data  of  the  type  described  earlier,  a  phenomenon  was  observed  that 
has  since  been  attributed  to  the  presence  of  inclusions  in  the  observed  liq¬ 
uid  droplets.  This  is  of  potential  interest,  because  droplets  with  internal 
bodies  may  occur  in  environments  such  as  hospital  buildings,  where  enti¬ 
ties  such  as  biological  organisms  can  remain  viable  over  long  periods  by 
encapsulation  within  liquids.  Similar  conditions  can  apply  in  the  dispersal 
of  biological  warfare  agents.  As  a  consequence,  there  is  interest  in  theoreti¬ 
cally  modeling  the  light-scattering  behavior  of  such  two-component  drop¬ 
lets,  with  the  objective  of  ultimately  establishing  a  means  of  recognizing 
their  presence  in  an  environment.  Reproducible  experimental  data  from 
droplets  with  inclusions  can  thus  provide  a  valuable  resource  for  verifying 
theoretical  inversion  methods. 

As  stated  earlier,  the  VOAG  droplets  are  normally  generated  into  a  dry  at¬ 
mosphere,  typically  <1  %RH.  When  20-|im  droplets  of  glycerin,  in  which 
water  is  soluble,  were  generated  into  a  moist  atmosphere  at  40  %RH,  no 
perturbations  of  the  expected  scattering  patterns  were  observed  (see 
fig.  8a).  However,  when  similar-sized  droplets  of  oleic  acid,  in  which  water 
is  insoluble,  were  generated  into  the  same  atmosphere,  light-scattering  pat¬ 
terns  of  the  type  shown  in  figure  8b  were  consistently  observed.  These  pre¬ 
liminary  scattering  data  were  recorded  at  a  flow  rate  of  1 1/ min. 

The  striking  feature  of  the  oleic  acid  pattern  is  the  displaced  secondary  scat¬ 
tering  fringes  superimposed  on  the  normal  concentric  rings  from  a  homo¬ 
geneous  sphere.  Figure  8c  shows  a  droplet  image  from  the  same  oleic  acid 
aerosol  (although,  for  the  reasons  given  in  sect.  2.2,  not  from  the  same  drop¬ 
let  as  shown  in  fig.  8b);  this  image  shows  the  presence  of  an  internal  dielec¬ 
tric  discontinuity.  The  positions  of  the  secondary  scattering  fringes  relative 
to  those  from  the  host  droplet  provide  good  evidence  to  support  the  hy¬ 
pothesis  of  an  inclusion  offset  vertically  upward  from  the  center  of  the  oleic 
acid  host,  and,  again,  this  is  supported  both  by  the  reproducibility  of  the 
images  of  the  droplets  (such  as  in  fig.  8c)  and  by  theoretical  analysis  [12]. 
This  upward  shift  of  the  inclusion  within  the  host  may  reasonably  be  as¬ 
sumed  to  be  due  to  the  density  difference  between  the  inclusion  and  host, 
as  is  the  case  with  water  and  oleic  acid  (0.998  compared  to  0.895  g/ cm3). 
The  higher  density  water  inclusion  moves  toward  the  top  of  the  host  oleic 
acid  droplet  when  the  droplet  is  accelerated  downward  within  the  airflow 
through  the  APS  delivery  nozzle. 
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Figure  8.  Scattering 
patterns  recorded 
from  droplets 
generated  in  40-%RH 
environment: 

(a)  20-jim-diameter 
glycerin  droplet  with 
no  inclusion  and 

(b)  20-pm  oleic  acid 
droplet  with  in¬ 
clusion.  Image 

(c)  shows  type  of 
droplet  giving  rise  to 
scattering  in  image 
(b). 
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These  results  are  preliminary,  and  the  exact  mechanism  by  which  the  water 
inclusion  is  created  is  not  yet  known.  One  hypothesis  is  that  water  vapor 
condenses  onto  the  original  droplets  and  dissolves  into  the  alcohol  present. 
Once  the  alcohol  has  evaporated,  the  water,  which  has  a  higher  surface  ten¬ 
sion  than  oleic  acid,  forms  a  single  droplet  enveloped  by  the  oleic  acid. 
Thermodynamic  modeling  of  such  a  system  is  complex,  but  has  been  the 
subject  of  recent  work  by  other  authors  [15,16].  For  smaller  host  droplet 
sizes,  the  secondary  scattering  becomes  distorted  to  yield  downward  bright 
arms,  as  shown  in  figure  9.  Theoretical  analysis  [12]  suggests  that  this  change 
in  scattering  is  commensurate  with  an  increase  in  the  size  of  the  inclusion 
relative  to  the  host. 


Figure  9.  Scattering 
pattern  and  image  of 
16-pm-diameter  oleic 
acid  droplet 
containing  inclusion 
with  greater  size 
(relative  to  host 
droplet)  than  that 
shown  in  figure  8. 
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6.  Discussion 


The  experimental  data  presented  in  this  volume  illustrate  the  complex  spa¬ 
tial  light-scattering  behavior  of  droplets  distorted  by  accelerating  flow  fields 
and  droplets  containing  inclusions.  These  data  are,  we  believe,  the  first  of 
their  type  to  be  reported.  The  reproducibility  of  the  experimental  data  makes 
them  a  valuable  resource  in  the  development  and  testing  of  theoretical 
inversion  models.  Such  models  may  ultimately  provide  a  route  to  the  rapid 
characterization  and/ or  identification  of  complex  droplet  morphologies, 
and  this,  in  turn,  could  lead  to  advances  in  instrumentation  for  droplet  aero¬ 
dynamic  sizing  and  for  the  detection  of  potential  bioaerosol  hazards  in  hos¬ 
pital  and  battlefield  situations.  In  volume  2  of  this  report,  numerically 
calculated  spatial  light-scattering  data  are  presented  that  indicate  that  theo¬ 
retical  inversion  of  the  light-scattering  data  is  achievable  for  a  variety  of 
droplet  morphologies  relevant  to  the  experimental  conditions  described  here. 
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